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This paper describes a simple deviee whieh permits mass eomparisons in air without appeal to the correction 
for air buoyaney. The deviee eonsists of a canister which is evacuated and weighed on a Laboratory balance with a 
mass inside. A second weighing of another mass in the evacuated canister provides the desired mass comparison. 
The method was used to determine the mass difference between two stainless steel weights of widely differing 
densities. With knowledge of this mass difference and of the volume difference one may, by a simple air weighing 
of the two objects, determine directly the density of the air in the balance case. Densities of air determined by this 
method were compared with those calculated from the barometric pressure, the temperature, and the relative 
humidity of the Laboratory air. The experimental and calculated values agree throughout to within 1.0 /xg cm -3 
(where the normal air density is about 1.2 mg cm -3 ). The calculated and experimental values of day-to-day 
fluctuations in air density agree to within 0.5 />ig cm -3 . 
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1 . Introduction 

In recent years, the authors have participated in two 
separate electrochemical studies of the Faraday constant. 
Koch and Diehl made a determination of the Faraday by the 
coulometric titration of a base, 4-aminopyridine [l]; 1 Bower 
and Davis have made a series of determinations by the 
anodic dissolution of high-purity silver into perchloric acid 
solution [2]. Both of these measurements were designed to 
provide values of the Faraday constant with uncertainties of 
about one ppm. A constant of such accuracy might make a 
contribution toward the resolution of the discrepancy between 
the values for the Faraday obtained directly (electrochemi- 
cally) and indirectly (by calculation from accurately known 
physical constants). For details on this last point the reader 
is referred to Taylor et al. [3], and Cohen and Taylor [4). 

The accuracy of an electrochemical measurement can be 
no greater than that of the associated mass measurements. 
The uncertainty in the definition and transfer of the electrical 
units into the Faraday laboratory in the Electrical Measure- 
ments and Standards Division, National Bureau of Standards 



1 Figures in brackets indicate literature references at the end of this paper. 



(NBS) is believed to be no greater than 0.3 ppm with respect 
to the NBS maintained units. To keep to 1-2 ppm the 
uncertainty of the measurement of the Faraday constant, the 
mass measurements, with all their corrections, must have an 
uncertainty of no greater than about 0.9 ppm, irrespective of 
the materials used. 

Of the various corrections to be applied to mass determi- 
nation, such as those for the sensitivity of the balance, the 
calibration of the standard weights, etc., the correction for 
the air buoyancy appears as the least amenable to direct 
determination. For the buoyancy correction, one almost 
always relies on a calculation of air density by some 
algorithm whose entries are such laboratory observables as 
the barometric pressure, the temperature, the relative humid- 
ity and, sometimes, the carbon dioxide content of the air. 

The general question of the limits of validity of the various 
closely related air density algorithms advanced in the scien- 
tific literature and the question of the direct precise deter- 
mination of air density in a balance chamber are matters of 
considerable moment at present to the metrological commu- 
nity. The urgency of these questions, among others, led the 
Bureau International des Poids et Mesures to hold a meeting 
on precision mass determinations in November 1976 at 
Sevres, France [5]. 
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Our study was undertaken with the aim of making direct, 
accurate measurement of the air density in order to confirm 
the air density algorithm which we used in computing the 
buoyancy corrections to the mass determinations required for 
the Faraday experiments. 

2. The Problem and the Method 

The evaluation of the buoyancy correction to weighings of 
4-aminopyridine, aqueous solutions and other low-density 
materials offers some difficulties if one uses high-density 
(platinum or stainless steel) standards. The density of 4- 
aminopyridine is about 1.27 g cm -3 . The weight standards 
are made of stainless steel, density about 8 g cm -3 . For an 
uncertainty of only 1 ppm in the mass of 4-aminopyridine or 
aqueous solutions, the density of air, about 0.0012 g cm -3 at 
room temperature, must be known with an uncertainty no 
greater than 0.08 percent. 

We have made a direct physical determination of the air 
density in a balance chamber and have compared the results 
with the version of the algorithm which is given below. The 
experimental method, based on Archimedes' principle, has 
been used successfully before, notably by Baxter [6] who 
used a globe to determine his buoyancy corrections, and by 
Jaquerod and Borel [7], who weighed a float in various 
samples of air under standard conditions in order to obtain 
the variation in density. To this end we have employed 
stainless steel artifacts to permit direct measurement of 
density of air in a balance. We have been anticipated in the 
use of this principle in recent years by Bowman et al. [8] and 
by Bruhlmans and Eschbach [9]. 

We have furthermore employed a method of determining 
the mass of these objects by comparisons in vacuo. Rather 
than using the expensive, but perhaps more accurate, method 
of placing a balance in an evacuable enclosure, we have 
constructed an evacuable enclosure which may be weighed 
on a chemist's analytical balance. 

Let M H and M s be the nearly equal masses of a hollow, 
low density weight and a solid, high density weight, respec- 
tively. Let the respective exterior volumes be V H and V s . Let 
M H and M s be weighed, successively, within an evacuated 
enclosure of mass M E and exterior volume V E . Calling the air 
density p, we have the force equation: 

g(M E + M H - pV E ) = (M E + M S - pV E + AM)g 

where AM represents the small difference indicated on the 
optical scale of the balance, calibrated with platinum 
weights, and g represents the local acceleration of gravity. 
Providing V E and p are constant throughout the procedure, 
the mass difference, AM = M H — M s has been obtained 
without appeal to an air buoyancy calculation. 

Once thus characterized, the two masses M H and M s , if 
their mean densities are sufficiently different, may be used 



to determine the density of the air in a balance case. If m is 
the difference in air weight (divided by g) between M H and 



Mc 



id if AV = V h 



m = (M s - pV s ) - (M H ~ pV H ), or (1) 

p = (AM + m)/AV. 

The values of p so obtained may be compared to those 
calculated from the atmospheric observables. 

The vacuum chamber is a stainless steel tube (fig. 1) with 
one end capped. To the other end is attached a stainless 
steel flange with a groove in which is set an 0-ring whose 
profile extends slightly above the groove. A separate flat 
stainless-steel cover with a stainless-steel vacuum valve 
welded to it reposes on the O-ring. A circular land on the 
underside keeps the cover centered on the cylindrical cham- 
ber. To the vacuum valve is welded a smooth stainless steel 
tube which serves to attach the whole chamber to a port on a 
pumping station. Removal of excess material from the valve 
enabled us to keep the mass of the whole assembly, with a 
17-g weight inside, within the 100-g capacity of our balance. 

The principal dimensions of the apparatus are approxi- 
mately as follows: length, bottom to flange, 85 mm; inside 
diameter of chamber, 27 mm; outside diameter of chamber, 
29 mm; outside diameter of flange and lid, 33 mm; overall 
height 127 mm. 

Two weights (fig. 1) were constructed of stainless steel, 
one hollow, the other solid. The hollow weight was con- 
structed by capping stainless steel tubing, 0.25 mm thick 
and 25 mm in diameter with stainless steel sheet 1 mm thick. 
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FIGURE 1. Vacuum weighing device consisting of a stainless steel canister 
and cover with vacuum valve. 

The low-density and high-density weights for determining density of air in a balance 
chamber appear at right. 
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The joinings were electron-beam welded and machined 
smooth, and the object was polished to a specular finish. The 
result was a weight in the form of a right circular cylinder 50 
mm long, 25 mm in diameter and weighing about 17 g (for 
ultimate use in a 20-gram capacity microbalance). The 
cylinder has an overall density of about 0.5 g cm -3 . The 
solid weight, a counter-weight to the hollow weight, was 
constructed from stainless steel bar stock 12 mm in diameter 
and is about 18 mm long. This weight was also given a mirror 
finish and during the polishing its mass was adjusted to be 
close to that of the hollow weight. The density of the solid 
weight is about 8 g cm -3 . 

Mass comparisons with the vacuum chamber were made 
on a commercial single-pan balance with air damping and 
optical scale. This balance had a capacity of 100 g and a 
resolution of 10 /xg. The air weighings of the two artifacts for 
checking the adequacy of the algorithm were all made on 
single-pan instruments with optical scale. These instruments 
all had a capacity of 20 g, and a resolution of 1 /xg. 

Humidity was measured by a simple psychrometer consist- 
ing of a wet bulb and a dry bulb thermometer, readable to 
0.1 C, ventilated by an electronic chassis fan. From the 
readings of this apparatus converted through Ferrels' tables 
[11] one obtains relative humidity or its equivalent as a term 
in the air density algorithm. The relative humity so obtained 
was compared to that indicated by a Dunmore-type electric 
hygrometer, calibrated at NBS. The two methods yielded 
values which showed a standard deviation from one another 
of less than one relative humidity percent. Such an under- 
tainty at 50 percent relative humidity is equivalent to an 
uncertainty of about 0.02 percent of the air density or about 
0.2 fig cm -3 . 

Throughout the comparison of mass and volume of the two 
artifacts runs the assumption that the volume of the hollow 
cylinder is constant near atmospheric pressure and that the 
volume of the chamber is constant when evacuated. The 
hollow weight has a permanent vacuum inside and the 
chamber is only weighed when evacuated. Shell and plate 
theory calculations [12] for the hollow cylinder indicate at 
worst a 30 ppm volume change as the pressure outside the 
cylinder varies from zero pressure to about one standard 
atmosphere. This is equivalent to about a 3 ppm variation in 
volume over the pressure range normally experienced in the 
Faraday laboratory at NBS. For the vacuum chamber, which 
is made somewhat more robust, these considerations hold a 
, fortiori. 

The comparison of the masses of the two objects may be 
briefly described. One of the objects is placed in the vacuum 
chamber and the chamber is pumped down to an internal 
pressure of 15 Pa (~100 mtorr) at which point the mass of 
air left in the chamber is about 10 fig, the resolution of our 
balance. The valve is closed and disconnected from the 
pump. The weight of the ensemble is compared against 



standard weights of stainless steel in the balance. All 
weighings are by double substitution with platinum sensitiv- 
ity weights. After the weighing, the chamber is connected to 
a source of dry nitrogen and is opened. The next object in 
the series is substituted for its mate in the chamber and the 
process is repeated. Weighings were carried out about 600 s 
after the evacuation in each instance. Therefore about 900 s 
elapsed between weighings. This period seemed to suffice for 
the dissipation of thermal gradients brought on by expansion 
of air through the valve hole and by handling. Further delay 
did not improve the precision. 

3. Results 



Table 1 shows a typical sequence of weighings of the 
evacuated chamber containing, successively, the hollow 
weight, the solid weight, and some platinum foil weights 



Tablk 1 

Successive weighings of evacuated chamber with artifact inside. Weighings 
were made at 15 minute intervals. 



Weighing Se- 
quence 


Platinum (g) 


Solid Stainless 
steel (g) 


Hollow Stainless 
steel (g) 


1 


94.56442 






2 




94.58228 




3 


94.56433 






4 






94.60929 


5 


94.56435 






6 




94.58233 




7 


94.56446 






8 






94.60937 


9 


94.56446 






10 




94.58236 




JJ 


94.56444 






12 






94.60942 


13 


94.56451 






14 




94.58246 




15 


94.56453 






16 






94.60943 


17 


94.56446 







which we had hoped would provide mass standards for use in 
vacuum. Figure 2 shows the data graphed. For weighings in 
vacuum the platinum results show high dispersion about a 
fitted line (cr,- = 48 fig), about 3 times the dispersion shown 
by the stainless steel weights. The slope in figure 2 arises 
from the slow rise of temperature in the balance case as the 
series of weighings progresses. The slope may be annulled, 
if desired, by treating the air as an ideal gas, and applying 
the correction for the exterior volume of the chamber. The 
apparent difference in slope is not statistically significant. 

The difference in mass of the two stainless steel artifacts 
was calculated in the following way. 
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Successive Weighings Of Evacuated 
Container With Artifact 







Hollow Stainless Steel Weight 


50 






94.60940 


_ 


• • 


grams 




• 


30 




• 

I I I 1 



50 

94.58240 _ 
grams 

30- 



4 8 12 16 

Solid Stainless Steel Weight 



I 



I 



50- 

94.56440 _• 
grams 

30 — 



6 10 

Platinum Weight 



14 



I I I I 



I I I I 



13 5 7 9 11 13 15 17 

Weighing Sequence 900 s Intervals 

Figure 2. Successive weighings of stainless steel and platinum weights in 
the evacuated canister. 



Straight lines were fitted to the two sets of data and <r pp , 
the standard deviation of a predicted point was calculated 
[10]. We chose the difference between predicted points at a 
point on the abscissa where the value of cr pp is nearly 
minimum for both sets. The minimum occurs near weighing 
number 9. We selected weighing number 9 and calculated 
the combined cr pp value of both sets. Thus 

average difference = 0.02699 g 
combined cr pp = 0.000017 g 

The uncertainty of the difference is thus only slightly larger 
than the 14 /xg precision of our 100 gram single-pan balance 
for a set of two weighings. We shall take 17 /xg as the 
uncertainty in the assignment of mass difference of the two 
artifacts. This quantity will be equivalent to 0.06 percent in 
the air density. 

The volumes of the two artifacts were determined by 
hydrostatic weighing in an inert fluorocarbon liquid with 
single-crystal silicon as a volume standard according to the 
elegant method of Bowman et al. [8]. 



4. The Air Density Algorithm - Comparison of 
Experimental and Calailated Values. 

A reluctance to accept without scrutiny any of the pub- 
lished air-density algorithms led us to reexamine closely the 
algorithm used throughout the work reported here and in 
references [1, 2]. We first assume that moist air is composed 
of two ideal gases: dry air and water vapor. Further, we 
assume that the mass of moist air in a given volume is simply 
the sum of the masses of dry air and water vapor as derived 
from ideal gas behavior. Finally, we account for the non- 
ideal behavior of moist air by introducing a dimensionless 
multiplicative function, the compressibility factor, which 
forces our results to agree with the work of Goff and Gratch 
[13], Hyland [14], and Hyland and Wexler [15] on the 
thermodynamic properties of moist air. The compressibility 
factor is nearly constant and differs from 1 by less than 0.05 
percent over the ranges of interest. 

Although the formulations differ, the assumptions of the 
algorithm we used are in fact those of the algorithm contained 
in the Smithsonian Meteorological Tables [16]. The Smith- 
sonian derivation lumps the difference in molecular weight 
between dry and wet air with the temperature, thereby 
creating a "virtual temperature." Following this compaction, 
the derivation again lumps the compressibility factor with the 
"virtual temperature" to create the "adjusted virtual temper- 
ature." 

In our formulation, the compressibility factors retain their 
identity. We calculated the density of air from: 



PoT„C da (0) 
C(T)PJ 



P - 



(M d 



M„ 



M A 



(2) 



vher 



P = p H ,{0) g(x) HO) 



(3) 



and the molar gas constant, R, is used in the calculation of 
Po, the density of dry air: 



Po 



M da P 
Caa(0)RT 



(4) 



In these equations the symbols have the following signifi- 
cance, and values [13], p is the density of dry air; T , the 
standard temperature, 273. 15K; C da (0), the compressibility 
factor for dry air at °C, 0.9994; P , the standard pressure, 
1013.25 millibars (101325 Pa); C(T), the compressibility 
factor of moist air at temperature F, 0.9996 at 20 °C and 
50% relative humidity; T is the absolute temperature, P is 
the pressure in millibars (10 2 Pa); M da is the effective 
molecular weight for dry air, 28.966 g mole -1 ; M w is the 
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molecular weight of water vapor, 18.016 g mole l ; and e, 
the saturation vapor pressure of water vapor in millibars, (10 2 
Pa). In (3), p Hg (0) is the density of mercury at °C, 13.5951 
g em -3 ; g(x) is the local acceleration of gravity on the second 
floor of the Metrology Building, NBS, 980. 102 cm sec -2 ; and 
l*(0) is the height in cm of the barometric mercury column 
corrected to °C. The value for the molar gas constant was 
taken from Cohen and Taylor [4]. 

The two masses were weighed on a 20-g single-pan 
balance (precision 1 fig). Values for air density in the 
balance chamber were computed by means of (1). These 
values were compared to air density values calculated from 
the air density algorithm using the barometric pressure and 
relative humidity of the laboratory atmosphere, and the 
temperature in the balance case. The comparison is listed in 
table 2. 

Also in table 2 are listed determinations of air density 
which were performed at the National Center for Atmospheric 
Research Boulder, Colorado, altitude above mean sea level, 
1864 m, and the Chemistry Department University of Wyo- 
ming, Laramie, altitude 2200 m. The measurements were 
carried out with balances and barometers similar to those in 
the Electrical Measurements and Standards Division, NBS, 
and with our own sensitivity weights and psychrometric 
apparatus. 

The readings taken at Boulder were marred by blustery 
weather. Also we think that the relative humidity reading 
may be too high for the conditions prevailing on the date of 
the measurement. The readings at Laramie were made in an 
underground laboratory. An unexpected shower occurred 
sometime during the measurements. However, since no 
abrupt changes in barometric pressure were observed, we 
have no reason to suspect the data. 



5. Discussion 

The data in table 2 indicate that once a value for AM is 
assigned by a vacuum mass comparison ( 1 ), the two artifacts 
will track the air density calculation to within about 0.05 
percent of the calculated value. The offset of —0.02 percent 
is very likely fortuitous, since the uncertainty of the assign- 
ment of AM, and hence of the offset, is ± 0.06 percent. 
Furthermore, the air density algorithm is as valid at high 
elevations as at NBS, Washington. There is, of course, no 
reason to expect a change with elevation, but the altitude of 
the two stations in Colorado and Wyoming is sufficient to 
show that uncertainties in Ap, table 2, due to differences in 
mean atmospheric pressure should be below 0.05 percent for 
the major metrological laboratories, most of which are low- 
altitude stations. 

In their paper on high precision density determination, 
Bowman et al. [8] report results on measurements of air 
density with hollow and solid weights. The random experi- 
mental uncertainty in the results was 0.03 percent to 0.05 
percent in the air density in the two separate series. Their 
measurements are independent of the air density algorithm 
since the value for AM was determined on a vacuum balance 
of commercial make | 18]. Bowman e/ al. |8] cite an average 
value for AM to the nearest microgram but give no error 
statement on this quantity. If we assume that their uncer- 
tainty was 50 /xg, the approximate resolution of the balance 
used [19], then their experimental uncertainty is about 0.14 
percent of the air density. 

Toropin and Snegov [20], in an original approach, adjust 
to identical weight in air an open bulb and a tube, both made 
of the same glass. The bulb is evacuated and sealed. The two 
devices now constitute a hollow and a solid weight of density 
0.125 g cm -3 and 2.0 g cm -3 respectively. Toropin and 



Table 2 
Comparison of the density of air, p, as measured with hollow and solid stainless steel weights, with that calculated from the algorithm. 



Hollow/solid 



Algorithm 



P 

(mg/cm 3 ) 



Diff. Ap(%) 



Bar. (mm Hg, °C) 



lC e(°C) 



t,«ry ( Q 



twet (°C) 



At NBS, Gaithersburg 



7 July 1976 


1.1608 


1.1603 


-0.04 


749.04 


24.8 


25.2 


18.0 


8 July 


1.1635 


1.1630 


-0.04 


750.34 


24.2 


25.3 


19.7 




1.1608 


1.1603 


-0.04 


749.44 


24.6 


25.3 


19.7 


9 July 


1.1683 


1.1676 


-0.06 


752.23 


24.2 


24.9 


18.0 


11 July 


1.1554 


1.1554 


0.00 


743.29 


23.6 


24.3 


18.4 


12 July 


1.1524 


1.1522 


-0.02 


740.20 


23.3 


23.9 


17.7 



At Boulder and Laramie 



10 Aug. 


0.9547 (B) 


0.9555 


+ 0.08 


615.11 


24.2 


23.9 


13.7 


11 Aug. 


0.9275 (L) 


0.9273 
Avg. difference 
Std. Dev. of Diff. 


-0.02 
-0.02 
±0.03 


592.30 


22.0 


21.5 


12.0 
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Snegov have therefore made the equivalent of a vacuum 
weighing. They found an offset of 0.01% in the density of air 
with respect to the calculation from their algorithm and a 
random scatter in p of about 0.04 percent, in agreement with 
the findings of Bowman et al. [8], and with our own. 

Jones [17] has recently reformulated an algorithm for 
calculating the density of air. He concludes from the 300 
ppm random uncertainty of his formulation and the 200 ppm 
systematic uncertainty that the uncertainty of mass transfer 
from a platinum-iridium to stainless steel kilogram is 30 /xg 
random uncertainty and 20 /xg systematic. Our values would 
yield by comparison a 70 tig overall uncertainty. 

Our experimental studies would assign a one standard 
deviation uncertainty of 0.08 percent to our calculation of 
the density of air. This amount would lead to uncertainty in 
assignment of mass values of only 0.08 mg per 100 cm 3 
difference in volume between objects. 

During the experiments we desired to alter the normal 
properties of the cylinders as little as possible when compar- 
ing them in vacuo. It was for this reason that we pumped no 
lower than 15 Pa. We presume that water absorbed on the 
stainless steel will be partly removed by the pumping, 
although a baking would be required to remove all of it, as is 
shown for stainless steel mass artifacts by Moret [19]. 
Yoshimori et al. [22] have shown that 0.2-0.7 tig of water 
per cm 2 may be removed from the surface of stainless steel 
weights merely by raising the temperature of the weights to 
100 °C. We must assume that we removed water randomly 
in our measurements in amounts comparable to those indi- 
cated by Yoshimuri et al. A difference of 24 cm 2 in the 
surfaces of our two weights could therefore yield an uncer- 
tainty due to absorbed water up to 17 /xg (approximately our 
experimental uncertainty). Kochsiek [23] has made a study 
of the effect of moisture upon the mass of foils of various 
metals used in the fabrication of mass standards. He finds 
that for mechanically or electrolytically polished stainless 
steel the amount of water adsorbed on the surface is on the 
order of 0. 1 /xg cm -2 at 50 percent relative humidity, a result 
which is consistent with Yoshimuri's. 

These results suggest that for more precise weighings, or 
for weighings at lower pressures than those described here, 
one would have to be very sure that the surfaces areas of the 
two weights were identical. 

These measurements were sufficient to assure the weigh- 
ings of any of our Faraday materials to within 0.8 ppm. We 
therefore carried the investigation no further. 

6. Conclusions. 

1. A convenient and inexpensive device can be con- 
structed for the direct comparison of two masses in vacuum. 
The device consists of a vacuum chamber which can be 
weighed on a balance. We believe that the uncertainty in our 
measurements (0.08% on the air density) can be reduced but 



physical and chemical evidence suggests that scrupulous 
care and perhaps chemical determination of surface moisture 
may be necessary for improvement of the results as much as 
tenfold. 

2. Our experiments and that of Toropin and Snegov lead 
us to conclude further that the various similar air density 
algorithms are valid to within 0.08 percent. This in turn 
leads to the conclusion that a mass determination of the 
following materials will have at most the relative uncertainty 
indicated (due solely to the air density algorithm) starting 
from platinum-iridium standards: silver, 5 X 10 -8 ; stainless 
steel, 8 X 10~ 8 ; water, 9 X 1(T 7 . 

It may be concluded that our weighings of 1 molar 
perchloric acid (density = 1.05 g cm -3 ) and 4-aminopyridine 
(density = 1.27 g cm -3 ) described above have an uncertainty 
of no more than 0.9 ppm from this cause. 

3. These measurements show that the assignment of the 
mass unit from one kilogram to another can be made with an 
uncertainty attributable to the air density algorithm of no 
greater than 0. 10 mg per 100 cm 3 in the volume difference of 
the kilograms. 
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